X=[.99*ones(1,6)];
Y=[-360 -340 -300 -200 -100 -50];
Z=[.84 .84 .84 .84 .84 .84];

X=[X .97*ones(1,6)];
Y=[Y -360 -340 -300 -200 -100 -50];
Z=[Z .84 .84 .84 .84 .84 .88];

X=[X .95*ones(1,6)];
Y=[Y -360 -340 -300 -200 -100 -50];
Z=[Z .86 .86 .86 .86 .889 .92];

X=[X .90*ones(1,6)];
Y=[Y -360 -340 -300 -200 -100 -50];
Z=[Z .9210 .9210 .9210 .9230 .947 .987];

X=[X .85*ones(1,6)];
Y=[Y -360 -340 -300 -200 -100 -50];
Z=[Z .9531 .9532 .9533 .9585 .9810 1.0150];

X=[X .8*ones(1,16)];
Y=[Y -360 -340 -320 -300 -280 -260 -240 -220 -200 -180 -160 -140 -120 -100 -80 -60];
Z=[Z .9788 .9788 .9789 .9792 .9797 .9805 .9814 .9825 .9840 .9865 .989 .9933 .998 1.0043 1.0131 1.0237];

X=[X .75*ones(1,12)];
Y=[Y -360 -350 -340 -320 -300 -280 -260 -240 -220 -200 -180 -160];
Z=[Z .9965 .9966 .9966 .9969 .9975 .9979 .9986 .9995 1.001 1.0022 1.0045 1.0069];

X=[X .70*ones(1,6)];
Y=[Y -360 -340 -300 -260 -200 -150];
Z=[Z 1.0087 1.0088 1.0093 1.0109 1.0142 1.0211];

X=[X .65*ones(1,5)];
Y=[Y -360 -340 -300 -260 -250];
Z=[Z 1.019 1.022 1.035 1.0696 1.0895];

clf
xlin=linspace(min(X),max(X),30);
ylin=linspace(min(Y),max(Y),30);
[XX,YY]=meshgrid(xlin,ylin);
ZZ=griddata(X,Y,Z,XX,YY,'cubic');
mesh(XX,YY,ZZ);
axis tight;hold on
plot3(X,Y,Z,'.','MarkerSize',15)
